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ABSTRACT 

We use hydrodynamic cosmologicai simulations to study damped Lya (DLA) and 
Lyman limit (LL) absorption at redshifts z = 2—4 in five variants of the cold dark matter 
scenario: COBE-normalized (CCDM), cluster- normalized (SCDM), and tilted (n = 0.8) 
VL m = 1 models; and open (OCDM) and flat (LCDM) Q m = 0.4 models. Our standard 
simulations resolve the formation of dense concentrations of neutral gas in halos with 
circular velocity v c > v c ^ res ~ 140 km s _1 for Q m = 1 and 90 km s" 1 for Q m = 0.4, at 
z = 2; an additional LCDM simulation resolves halos down to v Cjres ~ 50 km s _1 at 
z = 3. We find a clear relation between HI column density and projected distance to 
the center of the nearest galaxy, with DLA absorption usually confined to galactocentric 
radii less than 10 — 15 kpc and LL absorption arising out to projected separations of 
30 kpc or more. If we consider only absorption in the halos resolved by our standard 
simulations, then all five models fall short of reproducing the observed abundance of 
DLA and LL systems at these redshifts. To estimate the absorption from lower mass 
halos, we fit a power-law to the relation between absorption area a and halo circular 
velocity v c in our simulations and extrapolate using the Jenkins et al. (2001) halo mass 
function; we do not apply this method to the TCDM model because it has too few 
halos at the level resolved by our simulation. In the two LCDM simulations, for which 
DLA results agree well in the mass regime of overlap, the mean cross-section for DLA 
absorption is a ~ ir(0.3R v i r ) 2 , much larger than the simple estimate a ~ ir(0.1R v i r ) 2 
based on collapse of the baryons to a centrifugally supported disk (R v i r is the halo virial 
radius). The cross sections for LL absorption are a ~ ^(O.Gi?^) 2 , with a dependence 
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on numerical resolution at the ~ 25% level. Detailed examination provides further 
evidence of non-equilibrium effects on absorption cross-section: for example, individual 
absorbers can be slightly smaller in more massive halos because gas sinks deeper into the 
potential wells, but more massive halos nonetheless have larger average cross-sections 
because they are more likely to have multiple gas concentrations. Our extrapolation 
procedure implies that all four models are consistent with the observed abundance of 
DLA systems if the fitted a(v c ) extends to v c ~ 50 — 80 km s^ 1 , and they may produce 
too much absorption if the relation continues to v c < 40 km s~ x . Matching the observed 
abundance of LL systems requires absorption in halos down to v c ~ 30 — 50 km s -1 . 
Our results suggest that LL absorption is closely akin to DLA absorption, arising in less 
massive halos or at larger galactocentric radii but not caused by processes acting on a 
radically different mass scale. Robust tests of cosmological models against the observed 
amount of high column density absorption will require simulations of representative 
volumes that resolve halos at the low-mass limit where they cease to harbor high column 
density absorbers, 30 ^ v c & 60 km s -1 . 

Subject headings: quasars: absorption lines, galaxies: formation, large scale structure 
of the Universe 

1. Introduction 

Systems producing absorption in the spectra of distant quasars offer an excellent probe of 
the early Universe. At high redshifts, they easily outnumber other observed tracers of cosmic 
structure, including both normal and active galaxies. The interpretation of low column density 
quasar absorption systems has undergone somewhat of a revolution during the past several years, 
with the recognition that they may consist of gas aggregating into mildly nonlinear structures 
analogous in their dynamical structure to today's galaxy superclusters (Cen et al. 1994; Petitjean 
et al. 1995; Zhang et al. 1995, 1997; Hernquist et al. 1996; Miralda-Escude et al. 1996; Bi & Davidsen 
1997; Hui, Gnedin, & Zhang 1997). However, damped Lya (DLA) absorbers, with neutral hydrogen 
column densities Nm > 10 20 ' 3 cm -2 , are usually thought to be associated with the dense interstellar 
gas of high-redshift galaxies, based on several lines of circumstantial evidence: similarity between 
the column densities of damped systems and the column densities through typical spiral disks today, 
rough agreement between the total mass of atomic hydrogen in damped absorbers at z ~ 3 and the 
total mass of stars today (Wolfe & Prochaska 1998), measurements of radial extents > Wh^ 1 kpc 
in two DLA systems (Briggs et al. 1989; Wolfe et al. 1993), and direct imaging of a number of 
DLA hosts from ground based and HST observations (Rao & Turnshek 1998; Turnshek et al. 2000; 
Djorgovski et al. 1996; Fontana et al. 1996; Moller & Warren 1998; Le Brun et al. 1997). The 
nature of Lyman limit (LL) absorbers, with Nm > 10 17 - 2 cm~ 2 , is less well understood, though 
most models associate them with the outer regions of galaxies (e.g. Mo & Miralda-Escude (1996)). 



-3- 



Analytic studies based on the Press-Schechter (1974) formalism suggested that the abundance 
of DLA systems might be a strong test of cosmological models, potentially ruling out those models 
with little power on galaxy scales at z = 3 (Kauffmann & Chariot 1994; Mo & Miralda-Escude 
1994). The most sophisticated of these calculations, that of Kauffmann (1996), implied that the 
"standard" cold dark matter model (SCDM, with Q m = l,h = H /W0 km s _1 Mpc -1 = 0.5, and a 
power spectrum normalization as ~ 0.7) could account for the observed abundance of high-redshift 
DLA systems, with about 30% of the absorption at z = 2.5 occurring in galaxies with halo circular 
velocities v c > 100 km s -1 . Katz et al. (1996, hereafter KWHM) presented the first predictions 
of the amount of DLA and LL absorption based on 3-dimensional hydrodynamic simulations, con- 
cluding that these simulations of the SCDM model came within a factor of two of matching the 
observed DLA abundance but fell nearly an order of magnitude short of reproducing observed LL 
absorption. Ma et al. (1997) "calibrated" DLA estimates from collisionless N-body simulations 
against the KWHM SCDM simulations, then applied this calibration to N-body simulations of 
cold+hot dark matter (CHDM) models. They concluded that the CHDM scenario failed to repro- 
duce the observed DLA abundance even with a neutrino fraction as low as O v = 0.2, strengthening 
the earlier, analytic arguments, which focused on CHDM with Q. v = 0.3. 

Quinn, Katz, k Efstathiou (1996; QKE hereafter) and Thoul & Weinberg (1996) find that 
halos with circular velocities v c ^ 40 km s _1 are unlikely to harbor DLA absorbers. Gas in halos 
below this limit does not collapse sufficiently to shield itself from the UV background and reach 
the necessary HI column densities. The shortcoming of the KWHM calculation was that it could 
not include the contribution from DLA and LL systems below its resolution limit, corresponding 
to a halo circular velocity v c ~ 100 km s _1 . Consequently, the simulations themselves can only 
provide a lower limit to the total amount of DLA and LL absorption in the Universe. In Gardner 
et al. (1997a; GKHW hereafter), we addressed this shortcoming by combining the KWHM results 
with high resolution simulations of individual, low mass objects similar to those of QKE. We used 
these simulations to obtain a relation between absorption cross-section a and halo circular velocity 
v c , which we combined with the Press-Schechter halo abundance to compute the total DLA and 
LL absorption in the SCDM model. The correction for previously unresolved halos increased the 
predicted absorption by about a factor of two, bringing the predicted DLA abundance into good 
agreement with observations but leaving the predicted number of LL systems substantially below 
the observed number. In Gardner et al. (1997b; GKWH hereafter), we applied the a(v c ) relation 
derived for SCDM to other cosmological models, obtaining more general predictions for DLA and 
LL absorption under the assumption that the relation between halo v c and gas absorption cross- 
section was independent of cosmological parameters. 

In this paper, we present results of simulations of several variants of the inflation+CDM sce- 
nario (see, e.g. Katz, Hernquist & Weinberg 1999) and improve upon the GKWH results by us- 
ing these simulations to predict DLA and LL absorption in these models. We continue to use a 
Press-Schechter based extrapolation (with the mass function of Jenkins et al. 2001) to compute 
the contribution of smaller halos to DLA and LL statistics, employing an improved methodology 
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that significantly changes the GKHW predictions for absorption by low mass systems. Using an 
improved fitting procedure, we obtain more accurate error estimates of our fitted a(v c ) to the 
simulated data. We find that our largest error in estimating the universal amount of DLA and 
LL absorption arises from the uncertainty in the exact v c at which halos cease to harbor these 
absorbers. Given the large number of halos at v CjTn i n « 40 km s^ 1 , a small variation in the exact 
value or form of this cutoff leads to significant deviations in the estimation of total DLA and LL 
absorption cross sections. In light of these results, we find that we are not yet able to test the four 
cosmologies we consider against the observed DLA and LL abundances. Instead, we have adopted 
the approach of determining the value of v CjTn i n in each model that yields best agreement with the 
observations. 

The nature of the galaxies that host DLA systems has been a controversial topic for many 
years. Two competing hypotheses have defined the poles of the debate: the idea that most DLA 
systems are large, rotating gas disks (e.g. Schiano, Wolfe, & Chang 1990), and the idea that a large 
fraction of DLA absorption arises in dwarf galaxies (e.g. Tyson (1988)). The strongest empirical 
argument for the dwarf hypothesis is that some imaging studies reveal small galaxies near the line of 
sight but no clear candidates for large galaxies producing the absorption (e.g. Fontana et al. 1996; 
Le Brun et al. 1997; Moller & Warren 1998). The recent study of two DLA systems at z = 0.091 
and z = 0.221 by Rao & Turnshek (1998) and Turnshek et al. (2000) places especially stringent 
upper limits on the luminosities of the host galaxies. The strongest argument for the rotating disk 
hypothesis is the analysis of metal-line kinematics in DLA systems by Prochaska k, Wolfe (1997, 
1998), who consider a variety of simplified models for the velocity structure of the absorbers and find 
that only a population of cold, rotating disks with typical circular velocities v c > 200 km s _1 can 
account for the observed distribution of velocity spreads and for the high frequency of "lopsided" 
kinematic profiles. However, hierarchical models of galaxy formation predict that such massive disks 
should be rare at z ~ 3. In an important paper, Haehnelt, Steinmetz, h Rauch (1998) showed 
that hydrodynamic simulations of high-redshift galaxies could account for the lopsided kinematic 
profiles and large velocity spreads found by Prochaska & Wolfe (1997, 1998) even with halo circular 
velocities substantially below 200 km s _1 , because of large scale asymmetries and departures from 
dynamical equilibrium (see also Ledoux et al. (1998)). This result makes the ~ 100 km s _1 median 
halo circular velocities found by Kauffmann (1996) and GKHW for the SCDM model potentially 
compatible with the observed metal- line kinematics. Our present simulations do not yet have 
enough resolution for us to repeat the Haehnelt et al. (1998) analysis; we hope to do so with 
future simulations to carry out a statistical comparison between results from a randomly chosen 
cosmological volume and the Prochaska & Wolfe (1997, 1998) data. However, we can already extend 
the GKHW analysis to other cosmological models, predicting the fraction of DLA absorption arising 
in halos of different circular velocities. 

We also revisit an important issue explored by KWHM for the SCDM model, the predicted 
distribution of projected separations between DLA and LL systems and high-redshift galaxies. 
Given the number of recent attempts to directly image DLA host galaxies, our predictions will be 
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useful in testing the compatibility of the size and probable luminosity of our simulated DLA hosts 
with the imaging data. 

Section 2 describes the simulations and our analysis methods. Section 3 presents our analysis 
of the DLA and LL systems resolved by the simulations. Section 4 describes and applies our 
procedures for computing the contribution from unresolved halos. We discuss the implications of 
our results and present our conclusions in §5. 

2. Simulations and Methods 

2.1. The Simulations 

Our simulations follow the same general prescription as in GKHW, where a periodic cube whose 
edges measure ll.llft, -1 Mpc in comoving units is drawn randomly from a CDM universe and evolved 
to a redshift z = 2. We examine the effects of cosmology using five principal simulations detailed 
in the first five lines of Table 1, where as is the power spectrum normalization, h = km s _1 

Mpc -1 , Q m is the fraction of present-day closure density in matter, Qb is the fraction in baryons, 
and 17a = A/(3Hq), where A is the cosmological constant. The quantity n is the index of the 
inflationary fluctuation spectrum, with n = 1 corresponding to scale-invariant fluctuations. These 
are the same simulations presented by Katz et al. (1999), who studied the clustering properties of 
the galaxies. We will often refer to the three Q. m = 1 models (SCDM, CCDM, TCDM) collectively 
as the "critical" models and the two VL m = 0.4 models (OCDM, LCDM) as the "subcritical" models. 

The five principal simulations employ 64 3 gas and 64 3 dark matter particles, with a gravita- 
tional softening length of 5/i _1 comoving kpc (3/i -1 comoving kpc equivalent Plummer softening, 
l/i -1 physical kpc at z = 2). The particle masses are 1.5 x 1O 8 M and 2.8 x 1O 9 M for the gas 
and dark matter, respectively, in the critical models and 6.7 x 1O 7 M and 8.3 x 10 s M in the 
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Table 1: Model parameters. 
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subcritical models. These simulations were performed using TreeSPH (Hernquist k, Katz 1989), 
a code that unites smoothed particle hydrodynamics (SPH; Lucy (1977); Gingold & Monaghan 
(1977)) with the hierarchical tree method for computing gravitational forces (Barnes & Hut 1986; 
Hernquist 1987). 

The five principal simulations were done to study the effects of cosmology on DLA and LL 
systems. Because of uncertainties arising from resolution issues, we add to this study a further two 
"next-generation" simulations, L64 and L128, which were done in the same cosmology but with 
different mass resolutions to investigate the stability of our results with resolution. These were 
performed much more recently using PTreeSPH (Dave et al. 1999), a new parallelized version of 
TreeSPH, and their details are also given in Table 1. L64 is the same mass resolution as OCDM and 
LCDM, while L128 is a factor of eight greater in mass resolution (a factor of two greater spatially) 
and is valuable in examining absorbers in the lower-mass halos that the five principal simulations 
cannot resolve. 

Detailed descriptions of the simulation code and the radiation physics can be found in Hernquist 
& Katz (1989); Katz, Weinberg, & Hernquist (1996; hereafter KWH); and Dave, Dubinski, and 
Hernquist (1997). We only summarize the techniques here. For both simulation codes, dark 
matter, stars, and gas are all represented by particles; collisionless material is influenced only by 
gravity, while gas is subject to gravitational forces, pressure gradients, and shocks. We include 
the effects of radiative cooling, assuming primordial abundances, and Compton cooling. Ionization 
and heat input from a UV radiation background are incorporated in the simulation. We adopt 
the UV background spectrum of Haardt & Madau (1996), but reduce it in intensity by a factor 
of two at all redshifts so that the mean Lya forest flux decrement is close to the observed value 
given our assumed baryon density (Croft et al. 1997). We apply small further adjustments to the 
background intensity during the analysis stage to precisely match the Press, Rybicki, & Schneider 
(1993) measurements of the mean decrement (see Croft et al. 1997 for further discussion of this 
procedure). For example, at z = 3, the background intensity is reduced to 20% of the Haardt & 
Madau value, first by 50% during the simulation, then by a further 40% during post-processing. If 
we adopted a higher baryon density, f2& = 0.02/i" 2 instead of fi& = 0.0125/i~ 2 , then the background 
intensity matching the observed mean decrement would be a factor ~ 2.2 higher. We use a simple 
prescription to turn cold, dense gas into collisionless "star" particles. The prescription and its 
computational implementation are described in detail by KWH. Details of the numerical parameters 
can be found in Katz et al. (1999). 



2.2. Halo and Absorber Identification 

From the simulation outputs at z = 4, 3, and 2, we identify dark matter halos and the individual 
concentrations of cold, collapsed gas that they contain. We initially identify the halos by applying a 
friends-of-friends (FOF) algorithm to the combined distribution of dark matter and SPH particles, 
with a linking length equal to the mean interparticle separation on an isodensity contour of an 
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isothermal sphere with an enclosed average density contrast of 5 = 180. Then, the position of the 
most bound particle in each FOF-identified halo is passed on to the spherical density method (SO; 
Lacey & Cole (1994)) that calculates the sphere about the most bound particle that contains an 
overdensity of S = 180. The halos used in our analysis are those output from SO, and the circular 
velocities, which we denote as v c , are the actual circular velocities at the 5 = 180 radius (i?iso) of 
each halo. This method of characterizing halo mass has been shown to be the best for computing 
the halo mass function (Jenkins et al. 2001). 

To detect discrete regions of collapsed gas capable of producing Lyman limit and damped 
Lya absorption, we apply the algorithm of Stadel et al. (2001; see also KWH and http://www- 
hpcc.astro.washington.edu/tools/SKID) to the distribution of cold gas and star particles. SKID 
identifies gravitationally bound groups of particles that are associated with a common density 
maximum. Gas particles are only considered as potential members of a SKID group if they have 
a smoothed overdensity p g /p g — 1 > S v i r and temperature T < 30,000 K, and we discard groups 
with fewer than four members (we will apply a more stringent resolution cut below) . All of the gas 
concentrations found by this method reside within a larger friends of friends halo, even at z = 4. 
We match each absorber with its parent (SO) halo and discard halos that contain no absorbers. 
Including or excluding the "absorberless" halos in our mass function does not change the results 
above our resolution cutoff (explained below), since nearly all halos above our cutoff contain at 
least one absorber. 

We calculate the HI column densities for the halos by enclosing each halo within a sphere 
centered on the most tightly bound gas particle and of sufficient size to contain all the gas particles 
that might contribute to high column density absorption within the halo. We project the gas 
distribution within this sphere onto a uniform grid with a cell size of 5.43 comoving kpc, equal to the 
highest resolution achieved anywhere in the simulation, using the same spline kernel interpolation 
employed by the TreeSPH code for the hydrodynamics. For the L128 run, we use a pixel size of 2.715 
as the peak spatial resolution is a factor of two better in each dimension than the other simulations. 
Following KWHM, we calculate an initial HI column density for each grid point assuming that the 
gas is optically thin, then apply a self-shielding correction to yield a true HI column density (see 
KWHM for details). For each halo we compute the projected area over which it produces damped 
absorption, with Nm > 10 20 ' 3 cm~ 2 , and Lyman limit absorption, with iVHi > 10 17 ' 2 cm~ 2 . For 
simplicity, we project all halos from a single direction, although we obtain a similar fit of absorption 
area to circular velocity (see below) if we project randomly in the x, y, and z directions or average 
the projections in x, y, and z. Projecting a rectangular prism instead of a sphere yields the same 
results. To test for convergence, we reprojected several halos at 2 and 4 times smaller grid spacings 
and found that the cross section for DLA and LL absorption changed by less than 1% in the 
majority of cases and by at most 2.5%. 
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2.3. Numerical Resolution Considerations 

Our five principal simulations, which each contain 64 3 gas and 64 3 dark matter particles, 
lack the dynamic range needed to model simultaneously the full mass range of objects that can 
contribute to DLA and LL absorption. Simulations by QKE and GKHW show that halos with 
circular velocities as low as 35 km s _1 can host DLA absorbers, while photoionized gas is unable 
to collapse and cool in smaller halos. However, if we adopted a particle mass low enough to resolve 
35 km s _1 halos while retaining the same particle number, then our simulation volume would be 
too small to include a representative sample of more massive halos. 

In our analysis of the simulation results, we find that halos consisting of at least 60 dark matter 
particles nearly always (more than 98% of the time) contain a cold, dense gas concentration. Below 
this threshold, however, a substantial fraction of halos have no cold gas concentration. Furthermore, 
in our bootstrap analysis of the variance in the relation between halo v c and absorption cross-section, 
described in §4.1 below, we find much larger scatter about the mean relation for halos with fewer 
than 60 dark matter particles than for halos with more than 60 dark matter particles. 

To safeguard against additional systematic effects near the resolution boundary, we compare 
the high-resolution L128 run to L64. We find that the mass cut of M res = 60(m,dark + fnsPH) also 
allows the L64 halo properties to match smoothly with same mass halos in L128. Near perfect 
agreement for DLAs (cf. Figure 7 and later discussion) is found with a mass corresponding to 70 
particles. So to be as conservative as possible, we adopt M res = 70(mdark +rnspn) as an estimate 
of the limiting mass below which we cannot accurately compute the amount of absorption in a 
simulated halo. The 70 particle criterion is more conservative than the 34 particle criterion that we 
adopted in GKHW, and this change will affect our predictions for absorption in lower mass halos 
in §4 below. 

Applying this same mass cut to the L128 run, we find the limiting resolution to be roughly 
50 km s _1 at z = 3, which is unfortunately still above the mass at which halos cease to host DLA 
and LL absorption. Consequently, even the increased dynamic range of L128 does not allow us to 
model all absorbers. 

In the critical density models, the mass resolution limit is M res = 2.7 x 10 n M Q , corresponding 
to a circular velocity at the virial radius of v c ^ res = 140, 160, and 180 km s" 1 at z = 2, 3, and 4, 
respectively. In the subcritical models, M res = 8.2 x 10 10 M Q , corresponding to v c ^ res = 89, 100, 
and 112 km s _1 at z = 2, 3, 4. In L128, the mass resolution is a factor of eight better, hence 
M res = 1.0 x 1O 1O M and 

Veres — 50 km s ^ at z — 3. For our statistical analyses of the simulation 
results below, we always eliminate absorbers in halos whose total mass (dark matter plus baryons, 
with the spherical overdensity mass definition given above) is M < M res . Our quoted results apply 
only to halos above the resolution limit. In §4 we attempt to compute, as a function of v c , the 
contribution of halos with M < M res to the total amount of DLA and LL absorption by combining 
the Jenkins et al. (2001) mass function with our numerical results. 
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We assume throughout our subsequent analysis and discussion that our results for absorption 
in halos with M > M res are only minimally influenced by the residual effects of finite numerical 
resolution. We give evidence in §4 from the L128 run that the five principal runs are not influenced 
by resolution effects within a factor of 10 of their resolution cutoff. However, real absorption systems 
could have substructure that produces large fluctuations in HI column density on scales far below 
resolution limit of even our highest resolution simulation. In this scenario, the total amount of 
neutral gas in absorption systems would not be very different from our predictions, except to 
the extent that clumping shifts gas above or below the DLA/LL column density thresholds, but 
the distribution of column densities above the thresholds could be quite different. This issue will 
be difficult to address by direct numerical simulation alone because of the large range of scales 
involved. However, good agreement between predicted and observed column density distributions 
would support the contention that the absorbers do not have a great deal of substructure on scales 
below the simulation resolution limits. KWHM find good agreement between the predicted and 
observed shape of the column density distribution in the DLA regime, but a compelling case along 
these lines will require simulations that do resolve the full mass range of objects responsible for 
damped absorption, and which accurately resolve the low-end cutoff where halos cease to contain 
such absorption. 



3. Simulation Results 

3.1. Absorption in Collapsed Objects 

Figures 1 and 2 show the incidence of DLA and LL absorption in our five cosmological models: 
n(z) is the mean number of absorbers intercepted per unit redshift above the DLA (Fig. 1) or 
LL (Fig. 2) column density threshold. The numerical results for halos above the mass resolution 
limit M res are shown at z = 2, 3, and 4. Observational results for DLA absorption are taken from 
Storrie-Lombardi & Wolfe (2000; cf. also Storrie-Lombardi, Irwin, & McMahon 1996a; Wolfe et al. 
1995) and for LL absorption from Storrie-Lombardi et al. (1994). When comparing the subcritical 
models to the critical models, note that the subcritical models have lower M res and therefore sample 
the distribution of absorbers down to a lower mass cutoff, boosting the n(z) prediction relative to 
that of the critical models. Taken directly from the simulations and from halos only above M res , 
the values in these plots are hard lower limits to the predicted n(z). The limiting circular velocities 
v Cjres are below the value v c ~ 120 km s _1 inferred by Prochaska & Wolfe (1997, 1998) for typical 
DLA circular velocities based on a rotating disk model for metal-line kinematics, and even so the 
predicted number of DLA absorbers is usually a factor of two or more below the observed value. 
We conclude that if the inflationary CDM models considered here are even approximately correct, 
then the asymmetries and large velocity spreads found by Prochaska k, Wolfe must be a result of 
complex geometry and non-equilibrium dynamics, as proposed by Haehnelt et al. (1998). 

Figure 3 shows the distribution of impact parameters -Dproj, in physical units, between high 
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Fig. 1. — The incidence of DLA absorption in simulations of the five cosmological models listed 
in Table 1; n(z) is the number of systems with Nm > 10 20 ' 3 cm~ 2 intercepted per unit redshift. 
Simulation results are computed only for absorption in halos with mass above the mass resolution 
limits M res listed in Table 1. Observational data from Storrie-Lombardi & Wolfe (2000), shown 
with la (solid) and 2a (dashed) error bars, indicate DLA absorption by systems of all masses. The 
deficiency of absorption in the models may be explained partly or entirely by the contribution from 
halos below the simulations' resolution limits (see §4). 
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Fig. 2. — The incidence n{z) of LL absorption associated with halos above the simulation resolution 
limits. The format is similar to Fig. 1, but the threshold column density is now iVni > 10 17 ' 2 cm~ 2 . 
The upper error crosses represent the Lyman limit data of Storrie-Lombardi et al. (1994), with la 
and 2a abundance errors. The smooth curve shows their fitted power law. Observational values 
include systems of all masses. 
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Fig. 3. — The distribution of impact parameters -Dproj from high column density lines of sight to 
neighboring galaxies. The x-axis denotes the HI column density along the line of sight. The y-axis 
is the projected distance -Dproj (in physical kpc) to the galaxy in the simulation volume that lies 
closest (in projection) to the line of sight. The contours indicate the percentage of lines of sight 
containing an absorber of HI column density JVhi that have a galaxy whose center is within -D pro j . 
The thin contour levels are 10%, 25%, 75% and 90% with the dashed thick line denoting 50% and 
the solid thick line 99%. To construct the contours of the plot, lines of sight in the simulations 
were sorted into bins ANui = 0.5 dex and AD pro j = 2.5 kpc wide. 
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Fig. 4. — Average impact parameter (in physical kpc) vs. redshift. The bold lines are the average 
impact parameter (D pro j) of lines of sight with 10 17 < Nm < 10 17 ' 5 cm~ 2 . The lower thin set of 
lines is the same for 10 20 < Nm < 10 20 ' 5 cm -2 . The two curves at the top of the plot show the 
virial radius of a 150 km/s halo in a Vt m = 1 universe (solid) and £l m = 0.4 universe (dashed). 
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column density absorbers and the centers of neighboring galaxies. Specifically, the contour level 
represents the percentage of lines of sight of a given Ahi for which the closest simulated galaxy, in 
projection, lies within a projected distance -Dproj- Note that in this figure and in subsequent figures 
we represent distance in kpc and not /i" 1 kpc. Nearly all of the high column density systems in our 
simulations are associated with a galaxy, with the highest column density systems sampling the 
innermost regions of the galaxy and the lower column density systems occurring at larger impact 
parameters. At z = 2, nearly all DLA systems lie within 15-20 kpc of a galaxy center, and nearly 
all LL systems lie within 30 kpc. At higher redshifts, the most likely impact parameter increases, 
which could indicate a physical contraction of DLA systems as they age or could alternatively 
reflect the higher neutral fraction associated with a given overdensity at higher redshift (similar to 
the interpretation of evolution of the low column density forest given by Hernquist et al. (1996) 
and Dave et al. (1999)). This increase can easily be seen in Figure 4, which plots the mean impact 
parameter for systems at the DLA and LL cutoff and compares them with the virial radius of a 
v c = 150 km s _1 halo. The mass, and therefore size, of an isothermal sphere with a given circular 
velocity goes as (1 + z) -1 ' 5 . Not only does the mean absorption cross section increase at redshift 
z > 3, but the size of the parent halos decreases, meaning that at z = 4 the fraction of the area of 
the halo subtended by DLA and especially LL absorption is much larger than at z = 2. We will 
further examine absorber area vs. halo virial radius in Section 4.1. 

Figure 5 shows the fraction of critical density in cold collapsed gas, Sl ccg (solid line), and the 
fraction of the critical density in stars, VL± (dotted line), as a function of redshift in the various 
cosmological models. In subcritical models, we define Q x (z) = p x {z) x (1 + z)~ 3 /p c (z = 0), i.e., 
Q x represents the comoving density of component x relative to the critical density at z = 0. 
We obtain f2 ccs by integrating the column density distribution f(Nm) for all of the halos in the 
simulation. Error crosses show the values derived by Storrie-Lombardi, McMahon, & Irwin (1996b) 
from a sample of DLA systems. The largest observed column density for DLA systems in the 
Storrie-Lombardi et al. (1996ab) sample is 10 2L8 cm -2 , probably because higher column density 
systems are too rare to have been detected. For a more direct comparison to the data, we therefore 
compute an "observational" value, £l b s , for which we only count gas along lines of sight with 
-Am < 10 21 ' 8 cm -2 . The contribution to Q C cg from higher column density systems is generally 
small, but it is significant in the SCDM model. In all cases, we include only gas in halos with 

For all the models in Figure 5, 0* increases steadily as the Universe evolves. However, Q ccg 
remains constant to within a factor of two from redshift 4 down to redshift 2, indicating that 
additional gas cools and collapses to replace the gas that is turned into stars. Gas reaches higher 
densities in the the SCDM and CCDM models, leading to a larger difference between fi CC9 and f2 5 s . 
In nearly all cases, our "observed" cold gas densities £l b s fall at least a factor of two below the 
observational data of Storrie-Lombardi et al. (1996b), which could themselves be underestimates 
of the true cosmological values of Q ccg if dust extinction is important (Pei & Fall 1995). If any of 
these models are to be viable, a substantial fraction of the high redshift HI must reside in systems 



-15- 



below our resolution limit, an issue to which we turn in §4. 

To obtain the simulation values of Q ccg in Figure 5, we had to alter the procedure described 
in §2.2 for computing HI column densities. While our standard grid spacing of 5.43 comoving kpc 
is sufficient to resolve objects with HI columns of 10 20 ' 3 cm -2 and lower, a finer mesh is required 
to resolve the cold dense knots of gas at higher column densities, which contribute significantly 
to the Q ccg integral. As described in KWHM, we generate the initial HI map assuming complete 
transparency, then use the mass, HI mass fraction, and temperature of each grid cell to correct 
the HI column density for its ability to shield itself from the surrounding radiation. Typically, a 
high column density grid cell contains some regions where the hydrogen should be partly ionized 
and some where it should be completely neutral owing to self-shielding effects. At the standard 
resolution of 5.43 comoving kpc, our procedure may average the contributions of these two regions 
before the self-shielding correction is applied, resulting in a lower neutral column density than 
if the identical correction procedure were applied with a smaller grid spacing. This effect is not 
important for computing n(z), the number of systems with A%i above the DLA cutoff, but it can 
be important for computing the total mass density of neutral gas, £l ccg . We examined the effect by 
reprojecting some of the halos at z=2,3,4 in the SCDM model at 2 and 4 times the original spatial 
resolution. The original resolution of 5.43 kpc underpredicts the total HI in the simulation, while 
the cold collapsed gas mass in the 2X and 4X cases is nearly identical. Consequently, we regard 
the 2X case as numerically converged. 

Unfortunately, reprojecting all the simulation outputs at this higher resolution is not compu- 
tationally feasible. We therefore developed an approximate procedure based on the original grid 
spacing, calibrated against the few higher resolution SCDM projections. In grid cells where the 
self-shielding corrected HI column is greater than a threshold value Ahi, c , we treat as fully neu- 
tral all gas particles that contribute to that grid cell and meet the following criteria: temperature 
T < 30, 000K and gas density p g > ( 1000/1 77) p„j r (Clb/Q m ), where p v i r is the virialization overden- 
sity described in §2.2. For critical models, the density cut corresponds to 1000 In subcritical 
models, the density cut occurs at the same fraction of the critical density as in the Q m = 1 models. 
We find that for log A H i, c = (20.4, 20.7, 20.7) at z = (2,3,4) this procedure reproduces the SCDM 
high resolution values for O CC9 and Q b s to within 10%. 

3.2. Other Possible Sources of Absorption 

It is possible that some Lyman limit and/or DLA absorption originates from regions other 
than galactic halos. To investigate this alternative within our simulations, we project the entire 
simulation volume and compare the area of LL and DLA absorption to the sum of the absorption 
calculated by projecting each halo individually. In the analysis presented here, we use all halos 
that have at least one group identified by SKID as described in §2.2 (i.e. at least one concentration 
of cold gas that is gravitationally bound), whether or not the halo itself has M > M res . Above 
M = M res , 98% of the dark matter halos harbor at least one SKID-identified group. 
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Fig. 5. — Fraction of critical density in cold collapsed gas and stars in halos with v c > v c<res . The 
solid lines are the density parameter £l C cg in c °ld collapsed gas for each cosmological model. The 
dotted lines are f2*, the same quantity for the stars. The dot-dashed lines show "observational" 
values of 0, ccg , which include only lines of sight with HI columns iVm < 10 21 ' 8 cm -2 . The la error 
crosses are taken from Storrie-Lombardi et al. (1996b) and adjusted for the appropriate cosmology. 
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We have removed TCDM from the analysis in this section due to the extreme paucity of 
structure in the model. For the remaining four models, we calculate the total area subtended by 
DLA absorption in the halos with SKID-identified groups. Comparing this value to the total area 
subtended in an entire volume projection of each simulation at redshifts z = 2 and z = 4, we find 
agreement within 4.5% for all the models at both redshifts, and to better than 2% in five of the 
eight cases. We attribute the remaining differences to having more than one absorber along a given 
line of sight. Hence, all DLA absorption in the simulation occurs within halos with at least one 
concentration of cold, gravitationally-bound gas. 

In LL absorption, five of the eight outputs agree to better than 6% when compared in this 
manner. However, at z = 4 the results of volume projection and halo projection differed by 15%, 
30%, and 13% for SCDM, CCDM, and OCDM respectively. We took the worst case, z = 4 CCDM, 
and projected all the halos that contained at least 32 particles (gas + dark matter), whether or 
not they contained a SKID-identified gas concentration. When we sum the area subtended by LL 
absorption in these halos, we find that it now accounts for all but 1.2% of the LL absorption found 
by projecting the entire volume. Hence Lyman limit absorption still occurs exclusively in halos, 
but in this instance 30% of it occurs in halos in which our resolution of gas dynamics and cooling 
is only marginal and in which there are no SKID-identified gas concentrations. If our study had 
higher resolution, it is likely that some of these halos would have been able to form DLA systems as 
well. However, it is just to correct for these unresolved or under-resolved halos that we developed 
our Press-Schechter correction technique. The important conclusion is that all the Lyman limit 
absorption we find in the simulations resides within dark matter halos. If Lyman limit absorption 
were to occur outside galactic halos, it would have to be in regions that are much too small for us 
to resolve. 

To summarize, all DLA and LL absorption in our simulations occurs in dynamically bound 
dark matter halos, even below our resolution cutoff. At z = 2 and z = 4, in all four of the models 
tested, DLA absorption arises entirely in objects identified by SKID as bound concentrations of 
cold gas. LL absorption in the simulations also occurs exclusively in dark matter halos, although 
at z = 4 some of the gas within these halos is able to reach LL column densities without becoming 
a SKID-identified concentration. 

4. DLA and LL Absorption by Low Mass Halos 

4.1. Motivation from Higher Resolution Simulation 

We have so far focused on DLA and LL absorption in halos above our simulation mass resolution 
limits M res . However, if we want to test cosmological models against the observed incidence n(z), 
we must also consider the absorption that arises in lower mass halos, which are smaller in cross 
section but much more numerous. The L128 run is a factor of eight finer in mass resolution than 
the other simulations in this study, allowing us to examine trends in DLA and LL systems to lower 
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masses. Figure 6 compares, at z = 3, the circular velocity at R±sq of the halos in L128 with the 
cross section subtended by DLA absorption (left panel) and LL absorption (right panel) when the 
halo is projected. The number of vertices on each point indicates the number of SKID-identified 
concentrations of cold collapsed gas within the halo. We can see that the absorption characteristics 
of galactic halos are not well approximated by assuming a single galaxy per halo. Instead, the 
correlation of absorption cross section a(z,v c ) with halo mass seems to arise not from a single 
galaxy in each halo becoming larger as its parent halos increases in mass, but rather from more 
massive halos harboring more galaxies. Consequently, the multiple-absorber nature of halos is 
extremely important in modeling the connection between halo mass and absorption cross section. 
To approach the problem semi-analytically, it is necessary to model the full interaction history 
of the halos, as is done in Mailer et al. (2000). The absorption cross sections of the individual 
galaxies is virtually independent of halo mass. If any trend exists, it appears that a(z,v c ) in halos 
containing only one galaxy may actually decrease slightly in more massive halos. Higher mass halos 
have deeper potential wells, causing the concentrations of cold gas to contract more efficiently. This 
complex gas dynamical behavior demonstrates the value of a fully numerical treatment in modeling 
these objects. 

The dashed lines in Figure 6 show the area subtended by a face-on disk of radius R = 0.1R v i r 
where R v i r is the virial radius of an isothermal sphere with circular velocity v c at the virial radius. 
10% of the virial radius is the typical extent of a galaxy based on centrifugal arguments. The solid 
line in Figure 6 left (DLA) panel shows the area subtended by a face-on disk of radius 29%i?„; r , 
while in the right panel it denotes area a = ir (63% R v i r ) 2 ■ Although the solid lines were not fit to 
the data, one can see that the general trend of DLA and LL absorption is that the cross-section is 
roughly described as proportional to (30%i?t,j r ) 2 and (60%i?t,j r ) 2 respectively. 



4.2. Testing Extrapolations to Lower Mass Halos 

We would like to be able to extrapolate the contribution to the total incidence by halos with 
v c ,min < v c < v c . res . Given the results from Section 4.1, it is plausible to assume a power law fitting 
function api,(z, v c ) = Av^ . QKE find that a photoionizing background suppresses the collapse and 
cooling of gas in halos with circular velocities v c < v c m in — 

37 km s _1 . Thoul & Weinberg (1996) 
find a similar cutoff in simulations that are much higher resolution but assume spherical symmetry. 
As discussed in §2.1, v C:res is approximately 140 km s _1 in the critical models and 89 km s _1 in 
the subcritical models at z = 2 for the five principal runs. In L128, which has a resolution limit 
of 

Vc,res ~ 50 km s , we find no evidence of a photoionization cutoff; simulations with better 
mass resolution are required to detect it. For the a(v c ) dependence we find for resolved halos 
in our simulations, low mass halos dominate the total cross-section for DLA and LL absorption. 
Therefore, the predicted incidence n(z) depends sensitively on the assumed value of u c ,mm- 

Since we cannot robustly predict n(z) without exact knowledge of v C;m i n , we adopt the less 
ambitious goal of determining, for each cosmological model, what value of w c , m j„ yields a good match 
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to the observed values of n(z). Our approach to this calculation is similar to that of GKHW: we use 
our numerical simulations to calibrate a fit to the mean cross section for DLA (or LL) absorption 
of halos with circular velocity v c at redshift z, then integrate over an analytic halo mass function 
to compute n(z). 

Figure 7 compares a(v c ) in the L128 and L64 simulations at redshift z = 3. Figure 8 compares 
the cumulative incidence n(z,v c ), the total incidence from absorbers in halos at least at massive 
as v c , in L128 and L64. Note that in this Figure, the incidence is measured directly from the 
simulations and not by convolving an analytic mass function with o;pl(z,u c ) as is done later in 
this section. The runs were performed using the same cosmology but a factor of eight different 
mass resolution. In the absence of systematic resolution effects that may exist above our cutoff 
v c ^ res , the data from the L64 simulations should smoothly overlap with L128 above the L64 v C)res 
value. For the DLA case, although L64 appears to have more low-mass outliers than L128, the 
simulations agree quite well both in a(v c ) and incidence. In the LL case, however, we find that the 
we systematically underestimate the absorption cross section of LL absorbers in L64 by roughly 
25%. This leads to the cumulative incidence of the L64 simulation also being 25% less than the 
same v c in L128. Therefore, at the L64 resolution we are not resolving all LL absorption regions 
that exist inside halos with v c > v C)res . One possibility is that areas whose average HI column 
density is slightly below the LL cutoff (JVhi > 10 17 ' 2 cm~ 2 ) on ~ 5 kpc scales have smaller, patchy 
clumps with higher column density. Hence, the lower resolution simulations would systematically 
underestimate the LL absoption in these regions. For DLA absorption, however, we detect no 
signatures of numerical resolution artifacts. 

To find the best fitting function for absorption cross section, we bin halos in 0.05 dex increments 
in logw c , beginning with v c ^ res and subject to the constraint that there be at least 10 halos in 
each bin. We sometimes are forced to widen the bin size to satisfy the latter constraint. Let 
g dla = ot(y c ,z) denote the "cross section" of DLA absorption, i.e. the comoving area subtended 
by HI column densities 7V"hi > 10 20 ' 3 cm -2 when a halo is projected onto a plane. For the binned 
distribution of halos, we determine the log of the average halo DLA absorption cross section, 
log(o-DLA) for each bin. Then we calculate the statistical uncertainty of log(aoLA) in each bin by 
using the bootstrap method with 1000 random realizations of the data set of halos with M > M res . 
The distribution of halo cross sections (J dla at a given v c is approximately log- normal and hence 
best described by a Gaussian in log space. Since we will use the bootstrap errors in the next section 
to calculate confidence limits, which assume Gaussianity, we express the errors in log space. We fit 
the points log(aDLA) by linear least squares to determine the parameters A and B, the amplitude 
and index for the power law fitting function a-p-^{z,v c ) = Avf. The error crosses in Figure 7 show 
the mean absorption cross section in each bin of circular velocity at redshifts z = 3 for DLA and 
LL systems. The horizontal error bars show the width of each bin, and the vertical error bars show 
the la logarithmic uncertainty in {(Tola) for the bin determined by the bootstrap procedure. The 
solid error crosses are for L128 data and dashed error crosses for L64. The solid and dashed lines 
denotes the best fit to the L128 and L64 data respectively. For DLA systems, the fit to the L64 
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data is nearly identical to the L128 fit, showing that data at the resolution of the principal runs can 
be used reliably to estimate the absorption cross section below their resolution cutoff. The LL fit is 
systematically lower, reflecting the mismatch between L128 and L64 halo absorption. However, the 
two lines parallel each other meaning that the fitting procedure is robust and allows an accurate 
extrapolation to lower halo masses. 

It is important to note that we are not seeking to model the distribution of absorption cross 
sections in each bin, but only to characterize the mean cross section in each bin so that the total 
absorption can be accurately reconstructed from equation (1) below. The large spread and asym- 
metry in the distribution of absorption cross sections for individual halos is inconsequential for our 
purposes: api,(v Cl z) is the number which, when multiplied by the halo number density at v c , yields 
the total absorption at that v c that matches the total absorption present in the simulation. The 
bootstrap technique yields a robust estimate of the statistical uncertainty in this mean cross section 
caused by the finite number of halos in the simulation. Although the spread in absorption cross 
sections of individual halos may be large, their average absorption is a well determined quantity. 



Figure 9 shows the fitted relation a^{z, v c ) for each of 4 cosmologies and redshifts z = (2, 3, 4). 
Again we have removed TCDM from the analysis in this section due to the extreme paucity of 
structure in the model. Given so few halos above the mass limit, we felt it useless to attempt to 
fit to the data. The values of the fit parameters are detailed in Table 2. At z = 4, the fits for each 
model tend to follow the same slope (the exception being LL SCDM). At later redshifts, aph(z, v c ) 
evolves differently for different models. In general, LCDM tends to be among the steepest in all 
cosmologies. Interestingly, this steepness is not paralleled by OCDM at z = 2 where the cross- 
sections of the more massive halos in OCDM have decreased markedly. CCDM tends to be flatter 
than SCDM — the increased amplitude of structure apparently results in halos having smaller gas 
cross-sections. It is difficult to draw any general conclusions about the behavior of high column 
density absorbers in critical models vs. subcritical models, except that LCDM tends to have the 
highest cross-sections in massive halos. 

We are now in a position to correct the n{z) estimates from §3 to include the contribution from 
halos with M < M res . Our approach to this problem is similar to that of GKHW, though there are 
important differences of detail that make a significant difference to the end results, as we discuss 
later. We compute the number density of halos N(M, z) as a function of mass at specified redshift 
using the mass function of Jenkins et al. (2001; see their equation 9). Multiplying N(M, z) by our 
numerically calibrated functions aj>j,(v c , z), and integrating from v c to infinity, yields the number 
of DLA (LL) absorbers per unit redshift residing in halos of circular velocity greater than v c : 



4.3. 



Results 




(1) 
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where v' c is the circular velocity at Ri&o, the 5 = 180 radius, of a halo of mass M', and r is 
comoving distance (see GKHW for detailed discussion). If one takes the lower limit of the integral 
to be 

^c,mm? the minimum circular velocity for gas cooling and condensation, this yields tt-(z), the 
total incidence of DLA (or LL) absorption at redshift z. 

Figure 10 shows the results of this exercise. The curves show the cumulative incidence n(z, v c ) 
as a function of v c for redshifts z = (2,3,4). The error bars for SCDM and LCDM are shown 
at three representative locations and show the 1 — a error region resulting from the bootstrap 
uncertainty in the fits for api,(z,v c ). It comforting to note that the n(z,v c ) curve for LCDM, 
which is closest to the cosmology of L64 and L128, does indeed mirror the simulated n(z, v c ) shown 
in Figure 8. The cross-hatched region of Figure 10 denotes the 1 — a range allowed observationally 
by Storrie-Lombardi & Wolfe (2000; DLA) and Storrie-Lombardi et al. (1994; LL). To match DLA 
observations, contributions from halos with v c ^ 60 km s -1 are required. In the LL case, the 
minimum halo circular velocity is somewhat lower, more in the range of 40 km s _1 , 

although given the results from Figure 7, the estimate of n(z,v c ) may be depressed by ~ 25%. 
Raising n(z, v c ) by this amount would actually bring the minimum LL-harboring halo mass in line 
with the DLA estimate for most models. On the other hand, it is possible that LL absorbers could 
reside in halos of lower mass than DLAs. 

The results show that although halos in models like LCDM generally have higher absorption 
cross sections, the increased number density of halos in the critical models tends to give them more 
total absorption than the subcritical models. Although their a% values are close to that of SCDM 
(Table 1), and the growth factor reduction from z = 0toz = 2 — 4 is smaller in subcritical models, 
OCDM and LCDM also have redder power spectra and thus less power on these relatively small 
scales. The difference in power spectrum shape is especially important at the low-mass end of the 
mass function. 



5. Conclusions 

We can divide our conclusions into two classes, those that rely only on the results of our 
simulations, and those that rely on our extrapolation of these results via the Press-Schechter method 
(with Jenkins et al. 2001 mass function) to account for absorption by low mass halos. We will treat 
these two classes of conclusions in turn. 

5.1. Simulation Results 

Our five principal simulations resolve the formation of cold, dense gas concentrations in halos 
with v c > v c ^ res = 140, 160, 180 km s" 1 (89, 100, 110 km s" 1 ) at z = 2, 3, 4 in the critical density 
(subcritical) models. We employ a further 2 simulations done in identical cosmologies but with 
a factor of eight difference in mass resolution to examine resolution effects. The lower resolution 
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Fig. 6. — Comoving absorbing area (kpc 2 ) vs. circular velocity v c (km s _1 ) for halos in the LI 28 run 
at z = 3. The left and right panels show the area subtended by DLA and LL absorption respectively. 
The number of vertices in each data point corresponds to the number of gas concentrations in the 
halo, with the solid points representing halos containing a single absorber. In the left (DLA) panel, 
the solid line is the area subtended by a face-on disk of radius R = 0. 29 R v i r , where R v i r is the 
virial radius of an isothermal sphere with circular velocity v c at the virial radius. In the right (LL) 
panel, the solid line corresponds to a a face-on disk of radius R = O.Q3R v i r . The dashed line is a 
face-on disk of radius R = 0.1R v i r . 



Damped Lya 

z = 2 z = 3 z = 4 

Model B log A B log A B log A 

SCDM 1.196 0.722 0.904 1.777 1.049 1.605 

CCDM 0.936 1.232 0.601 2.275 1.036 1.594 

OCDM 0.398 2.403 1.207 1.067 1.238 1.254 

LCDM 1.247 0.764 1.569 0.325 1.151 1.508 

Lyman Limit 

z = 2 z = 3 z = 4 

Model B log A B log A B log A 

SCDM 1.079 1.585 1.443 1.144 2.213 -0.508 

CCDM 1.045 1.591 0.702 2.737 1.285 1.702 

OCDM 0.628 2.585 1.277 1.582 1.685 0.928 

LCDM 1.601 0.649 2.015 0.003 1.675 1.034 



Table 2: Fitted parameter values for ap^(v c , z) = Av, 
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Fig. 7. — Comoving absorbing area (kpc 2 ) vs. circular velocity v c (km s _1 ) for halos in the L128 
(crosses) and L64 (circles) runs at z = 3. As in Figure 6, the left panel shows DLA absorption 
and the right LL absorption. The horizontal components of the error crosses span the bins used 
to calculate bootstrap errors on the mean of a z (v c ), and vertical error bars show the 1-a limits. 
Details of the error estimation procedure are given in the text; note that the vertical error bars do 
not characterize the physical scatter in a(v c , z) but the statistical uncertainty in the mean value of 
a(v c , z) caused by the finite number of halos in the bin. The line is the best fit to the binned data. 
The solid line and error crosses denote the L128 data; the dashed line and error crosses show L64. 
The dotted line on the right panel corresponds to the best fit line for DLA systems (i.e. the solid 
line on the right panel). 
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run, L64, is the same resolution as the five principal simulations and has 

Vc,res equivalent to the 

subcritical (OCDM and LCDM) models. The higher resolution run, L128, h as v c res — 50 km s 
Our clearest conclusion is that absorption in halos above the circular velocity thresholds of the 
64 3 simulations cannot account for the observed incidence n(z) of DLA or LL absorption or for 
the amount of cold, collapsed gas, fl ccg , in observed DLA systems, for any of our five cosmological 
models (Figures 1, 2, 5). Higher resolution simulations are unlikely to change this conclusion, since 
clumping of the gas on scales below our gravitational softening length would tend to reduce the 
absorption cross section rather than increase it, unless this small scale clumping could produce 
neutral condensations in the outskirts of halos where we predict the gas to be mostly ionized. 
The evidence from L128 indicates that there are no resolution effects above v Ctres that affect DLA 
systems, although estimates of LL incidence in halos v c > v c/res may be underestimated by 25% in 
the principal simulations. 

Our models assume $7b = 0.0125/i~ 2 , and a higher baryon abundance (e.g. Buries & Tytler 
1998ab) might increase the predicted absorption. We have investigated SCDM models with different 
£lb values and find that higher Qf, leads to more absorption per halo as expected, but even a model 
with Qf, = 0.03125/i -2 = 0.125 has too little absorption at this circular velocity threshold to match 
the observations. We will report further results from this study in a future paper. 

If any of these cosmological models is correct, then a substantial fraction of high redshift DLA 
absorption must arise in halos with v c < 100 — 150 km s _1 . This conclusion appears consistent 
with the imaging of DLA fields, which often reveals no large, bright galaxies near the line of sight 
(Fontana et al. 1996; Le Brun et al. 1997; Moller & Warren 1998; Rao & Turnshek 1998; Turnshek 
et al. 2000). However, it implies that the asymmetric metal-line profiles found by Prochaska & 
Wolfe (1997, 1998) must be interpreted as a signature of non-equilibrium dynamics (Haehnelt et 
al. 1998) rather than smooth rotation. 

We find a clear and intuitively sensible relationship between high HI column density absorption 
and the proximity to galaxies (Figure 3). Damped systems typically lie within 10-15 kpc of the 
center of a host galaxy at 2 < z < 4, while lower column densities near the Lyman limit regime 
typically occur farther from the host galaxy. All DLA and LL absorption in our simulation occurs 
within collapsed dark matter halos. If it were to occur outside halos in the actual Universe, it 
would have to be on size scales smaller than we resolve. 

The stellar mass in our simulation is generally a steep function of time in the redshift range 
2 < z < 4, corresponding to a power law in z (Figure 5). The mass in cold collapsed gas, however, 
remains relatively fixed, indicating that the rate at which gas is converted into stars is roughly 
equal to the rate at which new gas cools out of ionized halos and condenses into galaxies. This 
result is expected if the star formation rate is an increasing function of gas density, as it is in our 
numerical formulation (KWH). 

In the first 3-d hydrodynamic study of high column density absorption, KWHM found that 
the predictions of n(z) from their simulations of the SCDM model fell a factor of two short of 
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the observed DLA abundance but a factor of ten short of the observed LL abundance. They 
speculated that the DLA shortfall could be made up by absorption in lower mass halos but that 
the LL shortfall might imply a distinct physical mechanism for the formation of LL systems, such 
as thermal instability on mass scales far below the simulation's resolution limits (Mo & Miralda- 
Escude 1996). It appears, however, that for the resolution of the KWHM runs and the principal 
simulations presented here, LL absorption at simulation scales may not have converged. If the 
resolution is increased by a factor of eight (as in the L128 run), LL absorption in halos over the 
same range in mass increases by 33%. This suggests the possibility that standard cosmological 
models can explain the observed LL systems with the physical processes that already occur in 
these simulations, albeit in halos somewhat below our current resolution limits. Thus LL and DLA 
absorption are closely related rather than physically distinct phenomena, with LL absorption arising 
preferentially at larger galactocentric distances and in less massive halos. We find no evidence in 
the simulations of LL absorption outside of galaxy dark matter halos. 



5.2. Absorption in Low Mass Halos 

In our simulations, the halo absorption cross sections a(v c , z) are determined by complex and 
competing physical processes. If we consider only halos that contain a single gas concentration, 
then the absorption cross section can actually decrease slightly with increasing circular velocity 
(solid points in Figure 7). However, more massive halos are more likely to contain multiple gas 
concentrations, with the net effect that a(v c , z) increases with increasing v c . At z = 2, the model 
with the weakest mass fluctuations (LCDM) tends to have high a(v c ) (Figure 10). This fact, 
and the trend for single-absorber halos, imply that DLA and LL absorption cross sections are 
substantially affected by non-equilibrium dynamics: absorbers get smaller if they have time to cool 
and condense in a quiescent dark matter potential well. Although others have argued that this 
behavior may be a numerical artifact (cf. Mailer et al. 2000), the agreement between our L64 and 
LI 28 runs and the appearance of the same trend in LI 28 suggests that it is not. Consequently, we 
suspect that this non-equilibrium behavior is a real feature of DLA and LL systems, possibly the 
geometric counterpart to the complex kinematic behavior found by Haehnelt et al. (1998). The 
physical complexity of a(v c , z) implies that an accurate fully analytic description of high column 
density absorption in CDM models will be difficult to achieve. Even the simple expectation that 
more small scale power produces a higher incidence of DLA and LL absorption does not always 
hold. 

In the L64 and L128 runs, for which DLA results agree well in the mass regime of overlap, 
the mean cross-section for DLA absorption is a « 7r(0.29i?„j,.) 2 , much larger than the simple 
estimate a ~ 7r(0.1i?„ r ) 2 based on collapse of the baryons to a centrifugally supported disk. For 
LL absorption, where we find that absorption in equal-mass halos is 25% lower in L64 than in 
L128 (which has a factor of eight finer mass resolution), the cross sections in L128 are described 
by a as 7r(0.63i?^ r ) 2 . 
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To estimate the amount of absorption in halos below the resolution limits of our simulations, 
we adopted a procedure similar to that of GKHW, using the numerical results to calibrate a(v c , z) 
and the Jenkins et al. (2001) mass function to compute the halo abundance. However, relative to 
GKHW we employed a much more conservative estimate of v C)res and an improved error estimation 
procedure based on bootstrap analysis instead of Poisson errors. These changes lead to superior 
a(v c , z) fits that generally increase the predicted amount of absorption in halos with v c <C v c ^ res . 
Our new results for n(z) in the SCDM model supersede those of GKHW, since our new procedures 
are certainly an improvement, and our results for n(z) in other models supersede those of GKWH, 
since in addition to these technical improvements we now have numerical simulations of these other 
models to constrain a(v c , z) for v c > 

Vc,res' The bootstrap procedure yields believable statistical 
uncertainties in the n(z, v c ) predictions. 

Taking our results and error estimates at face value, we find that four of the cosmological 
models that we consider are compatible with observational estimates of the incidence of DLA and 
LL absorption at z = 2, 3, and 4. What hinders us in better quantifying the total incidence in the 
Universe is our uncertainly in the estimate of the v cm i n , the circular velocity at which halos cease to 
harbor high column density systems. Previous studies (QKE; Thoul & Weinberg 1996) have found 
this cutoff to be approximately 40 km s _1 . However, due to the extreme number density of halos 
at this mass, a slight uncertainty in this cutoff results in huge uncertainties in estimating the total 
incidence in DLA and LL systems. Instead we determine, for each cosmology, the value of v Cjjn i n 
that best matches n(z) observations. Reproducing the data of Storrie-Lombardi & Wolfe (2000) and 
Storrie-Lombardi et al. (1994) requires u c ,mm ~ 60 km s _1 for DLA systems and v c ^ m i n ~ 40 km s _1 
for LL systems, with some dependence on cosmology and redshift (see Fig. 10). 

Since the DLA values of w c ,mm are above the expected threshold caused by photoionization, 
there is some risk that all of these models would predict too much DLA absorption in simulations 
that fully resolved the population of absorbing systems. A model with somewhat less small scale 
power, such as the lower amplitude LCDM model favored by recent Lya forest studies (McDonald 
et al. 2000; Croft et al. 2001), might fare better in this regard, perhaps matching the observed 
DLA abundance with a w c ,mm closer to the expected photoionization value. We are unable to make 
predictions for total absorption in the TCDM model with our current simulations because the 
paucity of structure above our resolution threshold makes our extrapolation procedure unreliable. 

Our current simulations provide a number of insights into the physics of DLA and LL ab- 
sorption in halos with v c > 100 km s _1 . Unfortunately, they also imply that robust numerical 
predictions of the incidence of high-redshift DLA and LL absorption will require simulations that 
resolve gas dynamics and cooling in halos with v c ~ 30 — 100 km s _1 , where our analytic modeling 
predicts a large fraction of the high column density absorption to occur. Simulations that resolve 
such halos exist (e.g. QKE; Navarro & Steinmetz 1997), but they do not yet model large enough 
volumes to predict statistical quantities like n(z). Achieving the necessary combination of resolu- 
tion and volume is challenging but within reach of current computational techniques. Simulations 
that meet these requirements will also teach us a great deal about the internal structure of more 
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massive DLA and LL systems and about the connection between these systems and the population 
of high redshift galaxies. 
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Fig. 8. — Cumulative incidence n(z, v c ) of DLA (left panels) and LL (right panels) systems for 
the raw simulated data from L128 (solid) and L64 (dashed) at z = 3. n(z,v c ) is the number 
of systems per unit redshift that are located within halos with circular velocity (at the 5 = 180 
radius) of at least v c . The cross-hatched region denotes the l-cr observed range of DLA absorption 
(Storrie-Lombardi & Wolfe 2000) and LL absorption (Storrie-Lombardi et al. 1994). 
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Fig. 9. — The best-fit absorption cross sections a(v c ,z) for DLA and LL systems in each model. 
The area is given in comoving kpc 2 . 
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Fig. 10. — Cumulative incidence n(z, v c ) of DLA (left panels) and LL (right panels) systems for 
each model at z = 4, 3, 2. n(z,v c ) is the number of systems per unit redshift that are located 
within halos with circular velocity (at the 5 = 180 radius) of at least v c . Representative error 
bars are plotted for the LCDM and SCDM models, with a horizontal offset of 0.01 dex applied to 
SCDM for clarity. The cross-hatched region denotes the 1-a observed range of DLA absorption 
(Storrie-Lombardi & Wolfe 2000) and LL absorption (Storrie-Lombardi et al. 1994). 



